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Abstract 

We present solutions to the Einstein-Maxwell system of equations in spherically sym- 
metric gravitational fields for static interior spacetimes with a specified form of the electric 
field intensity. The condition of pressure isotropy yields three category of solutions. The 
first category is expressible in terms of elementary functions and does not have an un- 
charged limit. The second category is given in terms of Bessel functions of half-integer 
order. These charged solutions satisfy a barotropic equation of state and contain Finch- 
Skea uncharged stars. The third category is obtained in terms of modified Bessel functions 
of half- integer order and does not have an uncharged limit. The physical features of the 
charged analogue of the Finch-Skea stars are studied in detail. In particular the condi- 
tion of causality is satisfied and the speed of sound does not exceed the speed of light. 
The physical analysis indicates that this analogue is a realistic model for static charged 
relativistic perfect fluid spheres. 



1 Introduction 

In the past many classes of exact solutions of the Einstein field equations have been found for 
static spherically symmetric gravitational fields with isotropic matter. A comprehensive list of 
Einstein solutions is provided by Delgaty and Lake [T] . These solutions may be used to model a 
neutral relativistic star as they are matchable to the Schwarzschild exterior at the boundary. In 
comparison fewer exact solutions of the Einstein-Maxwell field equations are known for static 
spherically symmetric gravitational fields with isotropic matter with nonzero electric fields. A 
recent review of Einstein-Maxwell solutions is given by Ivanov [2J. These solutions may be 
utilized to model a charged relativistic star as they match to the Reissner-Nordstrom exterior. 
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In modelling a charged relativistic sphere it is desirable to ensure that two general features 
are contained in the model. Firstly the model should be physically reasonable: the gravita- 
tional, electromagnetic and matter variables are continuous and well behaved in the interior, 
the interior metric matches smoothly to the exterior spacetime and causality is not violated. 
Secondly we should regain an uncharged solution (which should also satisfy the relevant physical 
conditions) of the the Einstein equations when the electrical field vanishes; a neutral star should 
be regainable as a stable equilibrium end state. The anisotropic charged model of Sharma et 
al |3j is an example that satisfies the stated criteria. An approach in satisfying the two criteria 
is to construct the model such that the limiting uncharged solution is a known exact solution. 
This is not easy to achieve in practice as the number of known exact uncharged solutions sat- 
isfying all conditions of physical acceptability are limited as established by Delgaty and Lake 

Our objective in the paper is to generate a new solution of the Einstein-Maxwell system that 
is physically acceptable and necessarily contains a neutral stellar model. The neutral stellar 
model is the Finch and Skea star [3] which satisfies all the requirement of physical acceptability. 
In §2 we express the Einstein equations for neutral matter and the Einstein-Maxwell system 
for charged matter as equivalent sets of differential equations utilising a transformation due to 
Durgapal and Bannerji |5]. We choose particular forms for one of the gravitational potentials 
and the the electric field intensity in §3. The condition of pressure isotropy becomes a second 
order linear equation in the remaining gravitational potential. We integrate the condition of 
pressure isotropy and consequently produce three classes of exact solutions to the Einstein- 
Maxwell field equations which can be written explicitly in terms of elementary functions as 
shown in §3.1, §3.2 and §3.3. We regain the uncharged exact solution found previously from 
the charged analogue of the Finch and Skea model in §3.2; we also demonstrate that the charged 
solution satisfies a barotropic equation of state. In §4 we comprehensively study the physical 
features of the charged Finch-Skea model. Graphs are generated for particular parameter 
values for the gravitational, electromagnetic and matter variables. We consider in particular 
the gravitational behaviour at the centre, and also produce a specific value for the the charge to 
radius ratio. We believe that this detailed analysis of the matter variables represents a realistic 
model of a compact charged object. Some brief concluding remarks are made in §5. 

2 Einstein- Maxwell equations 

We assume that the interior of a spherically symmetric relativistic star is described by the line 
element 

df ? = - e Mr) dt 2 + e 2X(r) dr 2 + ^2 + ^2 Q (1) 
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in Schwarzschild coordinates (t,r,8,(p) where the functions u(r) and A(r) are gravitational 
potentials. For neutral perfect fluids the Einstein field equations can be expressed as follows 
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^Kl-e- 2A )]' = P (2a) 



2u' 



-- (1 - e' 2X ) + — e~ 2A = p (2b) 

-^Ij/' + j/i + ^-j/X'--) = p (2c) 
r r J 

for the geometry described by (JIJ . We measure the energy density p and pressure p relative to 
the comoving fluid 4-velocity u a = e~ u 5^ and primes denote differentiation with respect to the 
radial coordinate r. We utilize units where the coupling constant = 1 and the speed of light 
c = 1. An equivalent form of the field equations is obtained if we introduce the transformation 

A 2 y 2 (x) = e 2v{r \ Z(x) = e" 2A(r) , x = Cr 2 (3) 

where the quantities A and C are arbitrary constants. Under the transformation (J3J) the system 
(0) has the form 

= L (4a) 

4Z^ + — = £ (4b) 
y x C 

AZx 2 y + 2Zx 2 y + (Zx - Z + l)y = (4c) 

where the dots denotes differentiation with respect to the variable x. 

A generalisation of the system (j3J) is the coupled Einstein-Maxwell field equations given by 

{ ■ E 2 x\ 
AZx 2 y + 2Zx 2 y + i Zx - Z + 1 - — J y = (5c) 

- = -{xE + E) (5d) 

where E is the electric field intensity and a is the proper charge density. When the electric 
field E = then the Einstein-Maxwell equations (0) reduce to the Einstein equations (@J) for 
neutral matter. The system of equations (jSJ) governs the behaviour of the gravitational field 
for a charged perfect fluid. The particular representation of the Einstein-Maxwell system as 
given in (jSJ) may be easier to integrate in certain situations as demonstrated by Thirukkanesh 
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and Maharaj JBJ. An interior stellar solution of (JHJ), for the line element (JIJ, should match to 
the exterior gravitational field described by the Reissner- Nordstrom line element. In terms of 
Schwarzschild coordinates the Reissner-Nordstrom solution has the form 

ds " = ~ (i - + 7?) dt " + ( x " V" + 1 dr2 + r '(^ 2 + sin2 ed ^ ( 6 ) 

where M and Q are associated with the mass and charge of the sphere respectively. 



3 Exact Solutions 

It is evident that the system of equations (JHJ) is underdetermined. Consequently we need 
to specify two of the variables in advance so that a solution can be obtained. A possible 
approach is to specify an equation of state relating p to p and additionally to specify one of 
the gravitational potentials. This leads to differential equations that are highly nonlinear and 
difficult to integrate. Another option, which is pursued here, is to postulate a form for the 
gravitational potential Z and to choose a form for the electrostatic field E. Our choice for the 
function Z(x) is contained in the form used by Thirukkanesh and Maharaj [0] and Maharaj 
and Mkhwanazi [ZJ which ensures that we regain as a special case models of fluid spheres with 
uncharged matter distributions analysed previously. 
The metric function Z(x) is chosen to be of the form 

Z = (l + x) n (7) 

where n ^ is a real number. (If n = then ()5a|) implies that p = —E 2 /2 which is negative.) 
The form (JJJ) ensures that the metric function e A behaves as 

e A = 1 + 0(r 2 ) 

near r = for suitable choices of n. In fact this is a sufficient condition for a static perfect fluid 
sphere to be regular at the centre as pointed out by Maartens and Maharaj jHj. We choose 
n = — 1 in (|7j) to complete the integration. Additionally we postulate the form 

= ax 

C (l + xf {) 

for E which depends on the real valued parameter a. The form (jHJ) is physically palatable since 
E 2 remains regular and positive throughout the sphere if a > 0. In addition the field intensity 
E becomes zero at the stellar centre and attains a maximum value of E = a/ aC/2 when 
r = 1/y/C. In what follows the Einstein-Maxwell equations (JHJ) are solved for all nonnegative 
values of the parameter a. 
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With the choices (J7J) (taking n = — 1) and (JBJ), equation (J5c|) reduces to 

4(1 + x)y - 2y + (1 - a) y = 



(9) 



It is convenient to categorise our solutions in terms of different values of a. We consider, in 
turn, the following three cases: 

a — 1, < a < 1, a > 1 
which generate classes of solutions to the Einstein-Maxwell system 



3.1 The case a = 1 

With the choice a = 1, equation 



assumes the simpler form 
4(1 + x)y - 2y = 



which is an ordinary differential equation of reducible order. It is easily integrated to yield 
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-ci(l + x) 3/2 + c 2 



where Ci and c 2 are constants of integration to be determined from physical considerations. The 
corresponding expressions for p and p can then be established via (}5aj) and (|5bj) respectively. 
The complete solution to the Einstein-Maxwell field equations © is then given by 



e 2A 

P_ 

C 

P_ 

c 

E 2 
a 2 



1 + Cr 2 

A 2 \ 2 Cl {l + Cr 2 f/ 2 + c 2 
Cr 2 + 6 



2(1 + Cr 2 ) 2 

ci(10 - Cr 2 )(l + Cr 2 f/ 2 - c 2 (2 + Cr 2 ) 
2(1 + Cr 2 ) [ci(l + Cr 2 ) 3 / 2 + 3c 2 ] 
C 2 r 2 



(1 + Cr 2 ) 2 
C 2 (3 + Cr 2 ) 



(10a) 
(10b) 

(10c) 

(lOd) 

(lOe) 

(lOf) 



(1 + CV 2 ) 5 

Note that this charged solution does not have an uncharged analogue as the electrostatic field 
intensity E cannot vanish (except at the centre). This effect essentially results from our condi- 
tion that a = 1. The line element for the solution (jlOjl is given by 

1 2 



ds 2 



-A 2 



^ Cl (l + Cr 2 ) 3 / 2 + c 2 



dt 2 + (1 + Cr 2 )dr 2 + r 2 (d6 2 + sin 2 



'IT 



In spite of the simplicity of this solution to the Einstein-Maxwell equations, it does not seem 
to have been published before. 
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3.2 The case < a < 1 



With < a < 1 equation is more difficult to solve. However it can be transformed to a 
standard Bessel equation. It is convenient to introduce the substitution X = 1 + x in to 
yield 



AX- 



d l Y dY 



-2— + (l-a )Y = 



(12) 



dX 2 dX 

where y(X) = Y. We now introduce a new function u(X) such that Y(X) = u(X)X m where 
m is a real number. Then the ordinary differential equation (|T2"|) becomes 



d 2 u du 
AX 2 —- + (8m - 2)X— + [Am 2 - 6m + (1 - a)X]u = 
dX dX 



(13) 



By introducing the further transform z = X 13 where (3 is a real number, (|13j) assumes the form 

A(3 2 X 2 ^ + [4/3(13 - l)X^ + (8m - 2)/3X p \ ^ + [4m 2 - 6m + (1 - a)X] u = 

This equation may be simplified by the choice [3 = \ and m — | which results in the differential 
equation 



, d 2 u du 



dz 2 ~^ dz ^ 



l-a)z z 

Now if we let (1 — a) x l 2 z = w then (|T4^) becomes 



u = 



(14) 



<i 2 w 



w 



dw 2 ^ dw ^ 



w 2 - 



u = 



(15) 



which is the Bessel equation of order |. 

The differential equation (|15J) has linearly independent solutions Js(w) and J_3(w) which 
are Bessel functions. The general solution to ()15|) can therefore be written as 



u = aJsiw) + bJ_3_{w) 



(16) 



where a and 6 are arbitrary constants. It is well known that the Bessel functions of half-integer 
order can be written in terms of the elementary trigonometric functions. In our case we obtain 
the following explicit forms for J a and J_3 in terms of elementary trigonometric sine and cosine 
functions: 



Ji(V^ — olz) 



sin \/T 



az 



J_3 — az) 



cos y/1 — az 



— cos Vl — 



az 



y/1 — airz 



\/l — az 



+ sin \/l — 
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Then the general solution to the field equation (J5cj) may be written as 



y{x) 



'1-a) 



-3/4 



[ci + c 2 Vl + xj sin - + x) 



+ c 2 



Ci \/i + xj cos a/(1 — a)(l + x) 



where we have introduced the new constants c\ = ay | and C2 = — oy | for simplicity. 

As is now determined the forms for p and p may then be established via (J5aj) and (|5b|) . 
The complete solution of the Einstein-Maxwell system (jSJ) for this configuration, in terms of 
the radial coordinate r, is thus given by the system 



e 
e v 

P_ 
C 
P_ 

c 



Vl + CV 2 
A 



Cl 



+ c 2 Vl + Cr 2 ) sin/(r) + (c 2 — c 1 vT+~CV 2 " ) cos/(r 



(1 - a) 3 / 4 L 
(2 - a)Cr 2 + 6 
2(1 + Cr 2 ) 2 

(/3(2 - a)(l + Cr 2 ) 3 / 2 + (a + 2 + 4a) (1 + Cr 2 ) + (3{a + 4d)Vl + Cr 2 - a) 
- ((a - 2)(1 + Cr 2 f /2 + /3(a + 2 + 4a) (1 + Cr 2 ) 



(17a) 
(17b) 

(17c) 



£7 



-(a + 4a) Vl + CV 2 - a/3) tan /(r 

-2(1 + Cr 2 ) 2 + Cr 2 - 1) - (J3 + Vl + Cr 2 ) tan/(r)) 

aC 2 r 2 



-1 



(1 + Cr 2 ) 2 
aC 2 (3 + Cr 2 ) 2 



(1 + Cr 2 ) 5 
where we have set 



(17d) 
(17e) 

(17f) 



[3 = a = vT 
c 2 



a 



1, /(r) = v/(l-a)(l + Cr 2 ). 



(c a + c 2 Vl + Cr 2 ) sin ^/(l - a)(l + Cr 2 ) + 



The line element for this solution has the form 
A 2 

~(1 - a) 3 / 2 

^c 2 - ciVl + Cr 2 ) cos v / (l-a)(l + Cr 2 ) 2 dt 2 + (1 + Cr 2 )aV 2 + r 2 (dfl 2 + sin 2 

We believe that this model is a new solution of the coupled Einstein-Maxwell system (J3J). It is 
easy to verify that on setting a = (i.e. a = 0) we regain the Finch and Skea (H solution for 
an uncharged sphere. From ()17ej) we have that E = as a = 0. The Finch and Skea solution 



7 



is then given by 

e A = Vl + Cr 2 



P 



P 



A (a + c 2 Vl + Cr 2 ) sin Vl + Cr 2 + ( c 2 - Cl Vl + Cr 2 ) cos Vl + Cr 2 
3 + CV 2 



;i + Cr 2 ) 2 

C (/Vl + Cr 2 + 1) + (/? - Vl + Cr 2 ) tan Vl + Cr 2 
~1 + CV 2 (/Vl + Cr 2 -l) - ((3 + Vl + Cr 2 ) tan VTTCV 2 



(19a) 
(19b) 

(19c) 
(19d) 



The uncharged solution (|19|) has been extensively studied by Finch and Skea and shown to 
be regular in the interior of the star and matches smoothly to the Schwarzschild exterior at 
the boundary. Furthermore, their solution has been shown to be consistent with the neutron 
star models proposed in the the theory of Walecka [5] . Our Einstein-Maxwell solution (JT7j) is 
a charged generalisation of the physically reasonable model (J19)) and does reduce to it when 
a = 0. 



3.3 The case a > 1 



We now briefly consider the case a > 1, so that 1 — a is negative, in the differential equation 
(151). It is convenient to introduce the new constant ip so that 



a 



As the integration of when a > 1 is similar to that in §3.2 we omit the details. The 
equivalent of equation (fHJ) is 

2" 



2 d u du 
dz 2 dz 



^ 2 z 2 + I - 



u = 



(20) 



If we introduce the new independent variable ipz = w then equation (1201) becomes 



dw 2 ^ dw 



d 2 u 



du 



w 2 + 



u = 



(21) 



We recognize ()21)1 as a special case of the modified Bessel equation of fractional order. It has 
linearly independent solutions, in terms of the modified Bessel functions 1 3 and J_s, which 

2 2 

may be expressed in terms of hyperbolic functions. These are given by 



lf)7TZ 



IpTIZ 



sinh^z) 

ipz 

smh{jpz) - 



+ cosh^z) 



cosh (ipz) 
ipz 
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The general solution to (}2"T]) is thus given by 

u = al-i{ipz) + bl_-i(ipz) 



(22) 



where a and b are constants. The equivalent of ()22j) is ()16|) for < a < 1. We can now obtain 
the gravitational potential y(x) in the form 

y{x) = ^CiVl + x — C2ipj sinh ypVl + xj + ^c 2 Vl + x — Ciipj cosh + xj 

which is the general solution of for a > 1. Note that we have introduced the new constants 
Ci = byf^ and c 2 = o<i/| in y(a;). 

The complete solution to the Einstein-Maxwell equations (0) can now be determined. In 
terms of the radial coordinate r it is given by 



Vl + Cr 2 



A 



ci 



Vl + Cr 2 - c 2 ^) sinh(^v / rTCV 2 ) 



+ (c 2 Vl + Cr 2 - dV>J cosh^Vl + Cr 2 ) 

(l-^ 2 )C 2 r 2 + 6C 
2(1 + Cr 2 ) 2 

(V^ 2 - l)C 2 r 2 - 2C 



(23a) 

(23b) 
(23c) 



2(1 + Cr 



2\2 



+ 



X 



1 + CV 2 
[ifj 2 (tanh(^Vl + Cr 2 ) + 0)] 



E 



a 



[((3ipVl + Cr 2 - 1) tanh(VVl + Cr 2 ) + (ipVl + Cr 2 - 0)] 
(^ 2 + l)C 2 r 2 

(1 + Cr 2 ) 2 
(ij 2 + l)C 2 (3 + Cr 2 ) 2 



(1 + CV 2 ) 5 

where we have put j3 = Ci/c 2 and a = ip 2 + 1. The line element for this solution is given by 



(23d) 

(23e) 
(23f) 



ds 2 



A 



CiVl + Cr 2 — C2ip ) sinh(ipVl + Cr 2 ) + ( c 2 Vl + Cr 2 — c\i\) ) cosh^Vl + Cr 2 



+ (1 + Cr 2 )dr 2 + r 2 {d6 2 + sin 2 



dt 2 
(24) 



It is interesting to observe that this solution is similar to our charged analogue of the Finch and 
Skea solution generated in §3.2. However in this case it is not possible to eliminate the electric 
field E (except at the centre), in order to obtain an uncharged counterpart, since ip 2 + 1 > 0. 
This means that this solution models a sphere that is always charged and hence cannot attain 
a neutral state. This feature is also shared by other charged solutions including that of Patel 
and Mehta [JU]. 
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In summary the Einstein-Maxwell models represented by the interior metrics (fTHj) and 
(|24j) are three new classes of exact solutions to the field equations which are all expressible in 
terms of elementary functions. 



4 Physical considerations: < a < 1 

We now study the physical properties of the solution (|T7|) corresponding to < a < 1 which is 
the uncharged analogue of the Finch and Skea neutron star model. Following Finch and Skea 
jl] we introduce the substitution 

v = Vl + Cr 2 

in our solution (j!7j) and obtain the following forms for the energy density p, the pressure p and 
the electric field intensity E: 

p (2-a)v 2 + a + 4 

C = 2^ (25a) 

V 



c 



[((2 - a)(3v 3 + (a - 2 + Ay/1 - a)v 2 + {a - 4 + 4y/l - a)(3v - a) + 

((a - 2)v 3 + (a - 2 + 4\/l - a)(3v 2 - (a - 4 + 4^1 - a)v - a/5) tan \A - av] 
x [-2v 4 ((pv -1)-(P + v) tan y/1 - av)] ^ (25b) 

in terms of the new variable v. 

The gravitational behaviour of this model is difficult to analyse because of the complexity 
of the expressions for p and p. In addition the solution is dependent on two parameters a and 
/?; it is necessary to select an appropriate value for one of the parameters and then determine 
the second parameter from physical considerations. Specifying a value for a is not a simple 
matter. Some values of a lead to unphysical behaviour (e.g. the value a — | results in the 
undesirable feature that po = for all choices of (3). Note that the subscript in p denotes the 
centre r = 0. In the present discussion we select the value 

17 

a = — 
81 

for the physical analysis of the matter and gravitational variables. This value of a chosen is 
consistent with the interval of validity for a (0 < a < 1) and it can be demonstrated that this 
particular value of a does not restrict the value of the parameter {3 in general. The matter 
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variables (|25|) assume the following simpler forms for this choice of a: 

p U5v 2 + 341 
C 



162v 4 



p (U5pv 3 + U3v 2 - 19(3v - 17) - (145i; 3 - U3f3v 2 - 19v + 17) tan (§v) 



C 
~C 



162t> 4 [(1 - /3v) + (J3 + v) tan (§«)] 



17(t; 2 - 1) 
81^ 



(26a) 
(26b) 



(26c) 



It is interesting to note that the Einstein-Maxwell solution (j2T?j) satisfies a barotropic equation 
V = p{p) which may be determined explicitly. The pressure p is given in terms of the energy 
density by 

^ = ( [l45p 3/2 + 143p - 19p 1/2 - 17] 

- [l45p 3/2 - 143p - 19p 1/2 + 17] tan fe 1/2 J ^ 



x 



( 162p 2 



(1 - pV2) + (1 + -1/2) tan ( ^1/2 



-1 



(27) 



where we have set 



145C + v/21025C 2 + 220968Cp' 
324p~ 



1/2 



for convenience. Thus (|2*7|). with a = 5?, represents a simple equation of state for a charged 



star. 



With the value of a = 5? the rate of change of the energy density, pressure and the electric 



field intensity are given by: 



dp 

~Cd~v 
dp 



-682 - 145w 2 



I- fiv + (P + v) tan ( -f 
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(28a) 



[435/3w 2 + 286w - 19(3 



-(435v 2 - 286/3u - 19) tan l-v J - -(U5v 3 - lA3pv 2 - I9v + 17(3) sec 2 ^ - ) v 

U5(3v 3 + U3v 2 - 19 Pv - 17 - (U5v 3 - U3pv 2 - 19v + 17/3) tan 

/g \ g /8^~'^^ 
4 - 5/3w + (4/3 + 5t>) tan ( -v ) + -v(p + v) sec 2 ( -v 



x 



1 - Pv + (p + v) tan ( -v 

9 



9 J 9 

1 2 N 



9 



1 dE y/a{2-v 2 



C dv v 3 \ / v 2 - 1 



(28b) 
(28c) 
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Using (J25|l we obtain the following expression for the speed of sound: 



dp 
dp 



1 - (3v + (/3 + v ) tan ( -« 

9 



[(435/V + 286w - 19/?) 



-(435w 2 - 286/3v - 19) tan i-v J - -(145t> 3 - U3f3v 2 - 19u + 17/?) sec 2 ^ ) v) 
145/3v 3 + 143v 2 - 19/3?; - 17 - (145v 3 - 143/?v 2 - 19v + 17/?) tan ( -u 
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4 - 5/3v + (4/? +5u) tan -u + -v(P + v) sec" | -r 



-290v 2 - 1364) 



(l-/3v + (fi + v) tan ( -u 



(29) 



For a regular model we require that ^ and ^ are well behaved in the interior of the charged 



(it! 



star. 



To ensure that the star is regular at r = it is important to examine the behaviour of our 
model at the centre of the star. The quantities (}26|) . (|28j) and (}29|) have the following forms at 
the centre of the star: 



Po 
C 

C 



dp 

dp 
Ccfa 
dp 
dp 



2.232/3-0.104 

J+ 9.620 
-827 



81 

363.79/3 2 + 1198.15/3 + 772.6 
(0.232/3 + 2.232) 2 
0.22/3 2 + 0.723/3 + 0.905 
(0.232/3 + 2.232) 2 



(30a) 
(30b) 

(30c) 

(30d) 

(30e) 



The quantity E vanishes at r = 0. In order to guarantee the positivity of the central pressure 
we obtain the following restriction on (3 from (j30b|) : 

P < -9.62 or /3 > 0.104 

The adiabatic sound speed condition < -j- < 1 yields the additional constraint: 

-4.06 < p < 5.9 

for the speed of sound not to exceed the speed of light. The intersection of the above inequalities 
is given by 

0.104 < p < 5.9 
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for a regular centre. It is clear that for this range of (3 the requirements po > 0, (^) < and 
(cfo)o < ^ are trivially satisfied. 

Now we are in a position to investigate the gravitational behaviour of our model in the 
interior of the star. The behaviour of the model is illustrated best in terms of graphs of the 
matter variables and the gravitational potentials. These graphs have been generated with the 
assistance of the software package Mathematica. Based on the above limits obtained on f3, we 
select the value 

(3 = 1. 

The graphs of the various quantities are on the interval 1 < v < 1.6. Figures 1-3 represent the 
behaviour of the energy density p, the pressure p and the electric field intensity E, respectively. 
In Figure 4 we have plotted * on the interval 1 < v < 1.6. The metric functions e 2v and e 2A 
are plotted on the same interval in Figures 5 and 6, respectively. 

4 

3.5 

3 

2.5 



o 

Q. 



2 

1.5 
1 

0.5 



1.1 1.2 



1.3 

v 



1.4 1.5 1.6 



Figure 1: A plot of the energy density p/C versus v = \ / l + Cr 2 . 




Figure 2: A plot of the pressure p/C versus v = Vl + Cr 2 . 



From the plots in Figures 1-2, we observe that the energy density p and the pressure p are 
positive and monotonically decreasing functions in the interior of the star. We observe from 
Figure 3 that the electric field is positive and monotonically increasing on 1 < v < \/2, has a 
maximum value at v = V2, and then decreases very slowly to the boundary. The potentials 
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Figure 3: A plot of the electric field intensity E/y/C versus v = y/1 + Cr 2 . 

1 



§ 0.6 
a, 

"a . 4 



0.2 



1 1.1 1.2 1.3 1.4 1.5 1.6 

v 

Figure 4: A plot of the speed of sound dp/ dp versus v = VI + Cr 2 . 



e 2u and e 2A are regular in the interior of the star as illustrated in Figures 5 and 6 respectively. 
Thus the quantities p, p, E, e 2v and e 2A are continuous, regular and well behaved throughout 
the interval 1 < v < 1.6. From Figure 2 we note that the pressure vanishes at approximately 



1.5252 



(31) 



which fixes the boundary of the charged star. The value (J31|) is consistent with the value of v 
for which p(v) = where p is given by (j30b|) . Of course the boundary of the star will change for 
other values of a and (3. A pleasing feature of this model is the behaviour of the speed of sound. 
It can be observed from Figure 4 that the speed of sound is always less than unity everywhere 
for 1 < v < 1.5252. Therefore the causality principle < dp/ dp < 1 holds throughout the star 
and the speed of sound is always less than the speed of light. Thus our solution ()17|). with 
a = |y an d /3 = 1, satisfies the requirements for a physically reasonable charged star. 

Utilising the information obtained from the graphs we are now able to examine the boundary 
conditions. In particular we have to match the interior ()18|) to the Reissner-Nordstrom exterior 
(jUj). With the help of v — 1.5252 = y/1 + CR 2 , the vanishing surface pressure condition 
p(R) = 0, yields the following approximate dependence of the constant C on the radius R of 
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Figure 5: A plot of the gravitational potential e 2y versus v = yl + C 
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Figure 6: A plot of the gravitational potential e versus v = VI + Cr 2 . 



the star: 



C 



1.33 



(32) 



The continuity of the electric field across the boundary condition E(R) = Q/R 2 gives the 
relationship of the total charge Q of the star, as measured by an observer at infinity, with the 
radius R for a = 17/81: 

VrfR 3 

^ 9{1 + R 2 ) 

Then we are in a position to find the charge to radius ratio 



Q 

-| = 0.26 
R 



(33) 



On matching of the interior gravitational potentials e 2v and e 2X from (fTHj) with the exterior 
Reissner-Nordstrom solution © we obtain 

1 



1 + CR 2 



729 
512 



Cl 



(1 + Vl + CR 2 ) sin (^VTTCR^j + 
(1 - Vl + CR 2 ) cos + OR 2 J 



2M Q 2 
~ ~R + R 2 



1 



(34a) 



1 + CR 2 



(34b) 
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On substituting (J32|) in ()34b|) we find the following values for the constants c\ and C2: 



ci = 0.05444 = c 2 
while (I34a|) gives another expression for C: 



2MR-Q 2 

R 2 (R 2 -2MR + Q 2 ) ^ 

in terms of Q. Note that on setting Q = we obtain the same value of C as in the uncharged 
model of Finch and Skea 0j. Therefore the introduction of the electromagnetic field has the 
effect of decreasing the value of C; consequently the radius R of the charged star is affected via 
the relationship (J32J) . Substituting for C from (J3~2j) and for Q from in we obtain the 
following mass-radius relationship: 

M , . 

— = 0.3 (36) 



Note that the Buchdahl limit 



M/R < 4/9 



for stability is applicable for neutral stars. We observe that this limit is not violated in our 
solution (|17|) in the presence of charge. In addition from (|33|) and (|36|) we obtain the following 
mass-charge ratio 

which satisfies the condition for equilibrium 

M 2 > Q 2 

of charged stars with nonzero pressure in general as pointed out by Cooperstock and de la Cruz 



5 Conclusion 

In this paper we obtained three categories of solutions for the Einstein-Maxwell system flH} 
corresponding to the electric field intensity (JSJ). In the first category of solutions E 7^ (except 
at the centre) and the solutions can be written in terms of elementary functions. In the second 
category of solutions it is possible for E = throughout the interior spacetime and we conse- 
quently regain the Finch and Skea solution [3] which is a good model for a neutron star. The 
charged analogue of the Finch-Skea solution is given in terms of Bessel functions of half-integer 
order. We demonstrated that this charged solution satisfies an explicit barotropic equation of 
state. In the third category of solutions E 7^ (except at the centre) and the solutions can be 
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written in terms of modified Bessel functions of half-integer order. A detailed analysis of the 
charged analogue of the Finch-Skea solution was carried out. For specific parameter values the 
energy density p, the pressure p, the electric field intensity E, the speed of sound dp/ dp, the 
gravitational potential e 2u and the metric function e 2A were plotted. The profiles of the matter, 
electromagnetic and gravitation variables suggest that these quantities are well behaved so that 
the charged analogue of the Finch-Skea solution describes a realistic charged stellar body. It was 
shown that the speed of sound is less than the speed of light and consequently causality is not 
violated. In particular we observe that the mass-radius ratio is within the Buchdahl limit for 
uncharged stars and the mass-charge ratio is consistent with the requirements of Cooperstock 
and de la Cruz [TT]. As this model contains the physically acceptable Finch and Skea solution 
and appears to be physically viable, a detailed analysis of the stability should be carried out; 
this is the subject of ongoing research. 
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